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Computer  Aided  Design  For  Reactive  Power  Compensation 

N.  Swain,  N.N.  Bengiamin  and  F.H.  Holcomb 


Abstract 

Presented  in  this  paper  is  a  computer  aided  design  tool  to  determine  the  size  and  location 
of  shunt  capacitors  for  the  purpose  of  reducing  power  losses  in  distribution  networks.  The 
developed  computer  program  (CapSize)  is  windows  based  with  several  design  options  that  permit 
practically  feasible  solutions.  The  implemented  algorithm  is  based  on  a  rigorous  optimization 
method  and  ac  power  flow  calculations.  The  algorithm  accounts  for  voltage  profile,  load 
variations,  load  transfer  switching,  observation  measurements,  and  limited  availability  of  data.  The 
results  of  a  test  case  are  presented  for  illustration. 

1.  Introduction 

The  usual  method  of  ohmic  loss  reduction  in  distribution  networks  is  through  the  use  of 
shunt  capacitors.  The  shunt  capacitors  reduce  the  reactive  component  of  power  flow,  and  hence 
lead  to  a  better  power  factor  and  voltage  profile.  Ohmic  losses  are  due  to  the  flow  of  real  and 
reactive  power  in  conductors.  The  former  can  not  be  reduced  since  real  power  is  the  useful  work 
producing  component  and  it  has  to  flow  the  length  of  the  feeder  from  the  energy  source  to  the 
■  load.  The  losses  due  to  reactive  power  flow  can  theoretically  be  reduced  to  zero  by  placing 
capacitors  at  each  node  with  the  size  of  capacitor  being  equal  to  the  reactive  component  of  load  at 
that  node.  Although  this  method  provides  the  minimum  reactive  system  losses,  it  is  not  practical 
because  the  cost  of  compensating  capacitors  may  exceed  the  benefits  gained  from  reducing  these 
losses.  This  may  also  lead  to  excessive  and  non  economical  capacitor  switching.  Therefore,  a 
compromise  solution,  that  will  address  these  issues  and  account  for  load  variations,  is  needed. 

Several  methods  of  loss  reduction  in  distribution  systems  are  reported  in  the  literature  [1, 
2].  These  techniques  vary  in  complexity  according  to  the  assumptions  made  and  system  details 
accounted  for.  Assuming  a  uniform  conductor  with  uniformly  distributed  loads  and  no  laterals,  a 
rule  of  thumb  (the  famous  2/3  rd  rule)  suitable  for  hand  calculations  was  possible.  With  the  advent 
of  computer  technology  this  limiting  assumption  became  unpractical  and  can  hardly  be  used  for 
any  realistic  situation.  The  more  advanced  tecAhiques  have  all  required  solving  a  non-linear 
optimization  problem  with  the  combined  demand  and  energy  costs  forming  the  objective  function. 

Heuristic  methods  are  often  suggested  to  reduce  computational  complexities.  Enumerative 
techniques  such  as  branch  and  bound  or  even  explicit  search  are  used.  More  rigorous  schemes  are 
of  the  iterative  type.  Researchers  did  not  spare  much  effort  formulating  and  solving  ohmic  loss 
reduction  problems  with  voltage  dependence  [3],  feeder  laterals,  switching  capacrtors  [4],  phase 
unbalance  [6],  voltage  regulators,  transfer  switches  automation  [7],  and  load  variations  included 
[1].  The  computation  time  escalates  rapidly  due  to  complexity  of  the  employed  search  schemes. 
The  authors  are  aware  of  only  one  commercially  available  computer  program  for  capacitors 
location  and  sizing;  it  is  of  the  explicit  search  type. 
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This  paper  presents  a  design  tool  that  permits  a  successive  evolution  in  design  to  satisfy 
ractical  constraints.  The  design  engineer  is  able  to  alter  the  ideal  theoretical  design  and  tailor  it 
according  to  his/her  particular  engineering  judgment  and  desired  outcomes. 


The  implemented  algorithm  is  an  extension  to  the  fundamental  work  by  Grainger  and  Lee 
n  6].  The  optimization  algorithm  utilizes  normalized  equivalent  feeder  data  to  easily  account  for 
non-uniformity  in  feeder  configuration,  and  ac  power  flow  to  account  for  voltage  variations.  The 
search  technique  is  iterative  and  it  progresses  to  satisfy  an  equal  area  criterion  for  current 
distribution  as  the  reactive  current  curve  shifts  with  the  addition  of  new  shunt  capacitors. 


2.  Design  Methodology 


The  original  work  by  Grainger  and  Lee  [3,  6]  assumed  radial  feeders.  This  assumption  is 
modified  in  this  work  to  include  laterals.  Each  lateral  is  treated  as  a  radial  feeder  for  which  the 
number  of  capacitors,  their  sizes  and  locations,  loss  reduction  and  savings  are  determined. 
Depending  on  the  savings  and  the  size  (a  user  defined  parameter),  either  all  the  loads  or  a  load 
equivalent  (difference  between  the  load  and  the  capacitor)  is  then  concentrated  at  the  lateral  feed 
point  on  the  main  feeder.  This  process  is  repeated  for  all  laterals  and  sub-laterals  until  the  main 
feeder  becomes  a  pure  radial  feeder.  This  equivalent  feeder  is  then  analyzed  for  reactive  power 
compensation.  The  algorithm  also  provides  the  user  with  the  option  to  determine  the  capacitor 
sizes  and  their  locations  for  an  equivalent  radial  feeder  where  all  lateral  loads  are  treated  as 
concentrated  loads  at  their  respective  lateral  starting  points.  This  is  useful  when  the  lateral  losses 
ture  very  small  and  an  approximate  practical  solution  is  desirable. 

Motivated  by  the  need  to  accommodate  load  transfer  without  voltage  profile  deterioration 
due  to  overcompensation  in  reduced  load  sections,  the  network  is  divided  into  zones.  The  zones 
are  determined  based  on  the  location  of  load  transfer  switches.  Assuming  open  transfer  switches, 
each  load  island  constitutes  a  zone.  Zone-1  is  the  zone  that  directly  connected  to  the  source. 


The  reactive  compensation  method  is  applied  to  all  zones,  and  is  earned  out  in  the 
following  four  phases: 

Phase-1  of  the  solution  solves  the  ac  power  flow  equations  to  determine  node  voltages,  branch 
losses,  branch  currents  and  real  and  reactive  power  at  the  substation. 

Phase-2  employs  the  normalized  feeder  and  optimization  principles  of  [6],  uses  the  substation 
power  and  reactive  current  calculated  in  phase- 1  and  accounts  for  load  variation  of  the  system  to 
determine  optimum  location  for  the  capacitors.  This  phase  may  produce  a  near  optimal  solution 
with  regard  to  capacitor  size. 

Phase-3  solves  ac  power  flow  equations  iteratively  to  determine  the  optimum  size  of  capacitors 
for  maximum  loss  reduction. 


Phase-4  determines  the  savings  due  to  power  and  energy  loss  reduction,  and  the  load  equivalent. 


104 


say  Qls,  (difference  between  total  reactive  component  of  the  load  and  the  total  capacitor  in  the 
zone).  Also,  the  total  real  component  of  the  load,  say  Pis,  the  zone  is  determined. 

For  the  closest  zone  to  the  substation  (zone-1),  consider  all  the  loads  originally  in  that 
zone  and  loads  Pis  and  Qls  from  the  other  zones  at  their  respective  branch  points.  Use  the  phases 
described  above  to  determine  the  capacitor  sizes,  location,  reduction  in  losses  and  savings. 

Since  each  zone  is  considered  individually  for  reactive  compensation,  the  shunt  capacitors 
will  naturally  be  located  within  the  zone  boundaries  allowing  for  transfer  of  compensation  with 
the  load.  Therefore,  an  acceptable  reactive  compensation  level  and  voltage  profile  are  always 
maintained. 


3.  CapSize 

Capsize  is  a  WINDOWS-based  computer  program  that  determines  the  optimal  number, 
size  and  location  of  shunt  capacitors  for  radial  distribution  feeders  without  and  with  laterals. 
Three  basic  design  options  are  available; 

Option-1)  Determines  the  required  number,  size,  and  locations  of  capacitors  yielding  the  most 
optimal  solution.  A  sub-optimal  design  with  one  capacitor  is  also  possible  using  this  option. 

Option72)  Determines  the  optimal  locations  for  capacitors  of  pre-determined  size.  This  option  is 
desirable  if  the  capacitor  sizes  obtained  from  option- 1  are  not  in  the  practical  range. 

Option-3)  Determines  the  optimum  size  of  capacitors  at  pre-detemnned  locations.  This  sub- 
optimal  design  is  useful  if  option-1  produced  an  excessive  number  of  capacitors  or  the  determined 
locations  are  undesirable. 

The  design  engineer  may  use  these  options  in  any  order  to  meet  the  particular  objectives 
of  the  system.  Besides  these  three  options  CapSize  has  also  the  following  features; 


•  graphs  the  current  distribution  of  the  feeder  with  and  without  laterals. 

•  displays  and  prints  a  brief  summary  and  a  detailed  report. 

•  edits  the  data  files. 

The  three  design  options  share  some  common  schemes  in  their  algorithm.  These  schemes 
will  be  presented  first  and  the  individual  algorithms  will  follow. 

Common  Schemes: 

a.)  AC  Power  Flow  (Subprogram-APF) 

1.  It  determines  voltage  at  each  node,  current  and  losses  in  each  branch  of  the  feeder  using  power 
flow  equations;  determines  the  total  real  power  loss  (say  Pmax)’  calculates  the  maximum  possible 
real  power  loss  reduction. 
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(1) 


Pe'^max  '  ^min 

where,  Pniin  is  the  ohmic  loss  due  to  the  flow  of  real  power  component, 
b)  Feeder  Parameter  Normalization  (Subprogram-FPN) 

1.  It  normalizes  the  feeder  parameters  and  calculates  the  normalized  current  as  follows[6]: 

i.  Choose  the  resistance  r*  ohm/unit  length  of  the  first  section  of  the  feeder  as  the  base  for 
normalizing  the  uniform  feeder's  impedance. 


ii.  Calculate  the  equivalent  length  "Lui"  for  each  section  (i)  using 

Ljj|~(Li/ri)/r*,  i— l,2..k, 

iii.  Calculate  the  total  length  of  the  equivalent  uniform  feeder  "Lu"  using 

k 

Lq  =  Z  Lui 

i=l 

iv.  Calculate  the  equivalent  normalized  length  "Lni"  for  each  section  using 

Lni  ~  I^ui^u 

y.  Calculate  the  resistance  per  unit  length  "mo"  for  the  normalized  feeder  using 
k 

rno  =  Z  Lj*ri 

i=l 


(2) 

(3) 

(4) 

(5) 


vi.  Calculate  the  normalized  current  in  each  section  of  the  feeder  by  dividing  the  actual  feeder 
current  in  each  section  (ficni)  by  the  base  reactive  current,  (Qo^^-^^o) 

vii.  Increase  the  reactive  component  of  the  load  in  each  section  by  the  corresponding  section 
reactive  losses  calculated  in  (a). 

2.  Create  the  data  file  with  the  normalized  parameters. 


3.  Calculate  a  as  follows: 

a  =  (Kp  +  Ke  *  T  *  Lf)/(Kp+Ke  *  T)  (6) 

This  parameter  accounts  for  the  load  factor  effect,  where  Kp=  demand  charge  in  $/kw/yr, 
Ke=energy  charge  in  mills/kwh,  T=time  in  hrs,  and  Lf=annual  reactive  load  factor. 
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Individual  Algorithms 


a)  Option-1:  Optimum  Number,  Location  and  Size 

Step  -1.  Call  Subprogram-APF  to  determine  the  voltage  at  each  node,  current  and  losses  in 
each  branch  of  the  feeder. 

Step-2.  Call  Subprogram-FPN  to  create  the  normalized  parameter  data  file. 

Step  -  3.  The  objective  of  the  capacitor  placement  is  to  maximize  the  net  dollar  savings  resulting 
from  both  power  and  energy  loss  reductions  and  is  given  by, 

n 

S=  KpLP  +  KgLE  -  KcSQci  U) 

i=l 

where  LP=total  real  power  loss  reduction,  LE=  total  energy  loss  reduction,  and  Qci=  capacitor 
size,  n  is  the  number  of  capacitors,  and  Kq  =  cost  of  capacitor  in  $/  kvar.  Assume  —  0  and 
invoke  the  Equal  Area  Criterion  [6]  of  step-4  to  determine  the  capacitor  locations  and  sizes  to 
minimize  power  and  energy  losses  and  maximize  S.  This  is  the  near  optimal  solution. 


Step  -  4:  Equal  Area  Criterion: 

^Determine  the  optimum  locations  "hj"  and  size  "Qcj"  where  hj  is  the  distance  of  capacitor  j 
'  from  the  source'  Capacitor  number  "1"  is  the  farthest  from  the  source.  This  is  done  as  follows. 

4a.  Select  an  arbitrary  value  for  hj;  e.g.  the  end  of  the  feeder.  This  value  determines  the  farthest 
permissible  location  from  the  source  where  capacitors  can  be  located.  Apply  the  folio-wing  steps 
for  i=l,2, . ,n. 

4b.  Calculate  Ici,  the  current  due  to  capacitor  "i"  using 
i-1 

Ici=2{aI,F(hi)-ZIcj}  <*) 

j=l 

where  F(hi)  is  the  normalized  reactive  current  at  hj.  The  "2"  multiplier  is  relevant  to  the  popular 
2/3  rule  in  reactive  power  compensation.  Thd  capacitor  size  Qd  is  calculated  using  ^3  Vhi  ^ci- 

4c.  Find  "gi"  such  that 
i 

F(giH  z  IcjV(«*Is) 

where  Ig  is  the  peak  reactive  current  injected  to  the  feeder  at  the  source,  and  F(gi)  is  the 
normalized  reactive  current  at  location  "gi"  on  the  equivalent  feeder. 

4d.  Shift  the  normalized  feeder  reactive  current  curve  down  such  that  F(gi)=0,  then  find  hi+i 
which  produces  an  area  under  the  curve  equal  to  area  "A”  above  the  curve  that  is  bounded 
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by  hj  and  gi.  Figure  1.  If  hj+j  is  a  viable  location  on  the  feeder,  increment  "i”  and  go  to  4b. 


F(x) 

1.0 


Legend: 


h2  I 
Area  B1 


Original  Normalized  Reactive  Current 
Shifted  Normalized  Reactive  Current 


gi  T 

Area  A1 


■X=l 


X 


Figure  1:  Illustration  of  Equal  Area  Criterion. 


4e.  If  hn+i=0  and  An=Bn,  then  stop,  a  solution  has  been  found.  Otherwise  ,  reset  "i"  to  1,  reduce 
hi  if  Bn>An  and  increase  h]  if  B^  <  Aji,  then  return  to  4b.  It  must  be  noted  that  if  hj  is  a  step 
feed  point,  the  algorithm  has  the  flexibility  of  adjusting  F(hi)  to  satisfy  the  equal  area  criterion. 

Step-5.  Convert  the  normalized  locations  hj,  i=l,2,...,n  to  their  corresponding  actual  locations 
along  the  feeder. 

§|tep  -  6.  Call  subprogram  APF.  Calculate  the  new  losses,  voltages  and  reactive  currents.  Go  to 
step  8. 

Step  -  7.  Increase  the  capacitor  sizes  by  equal  amount  and  run  Subprogram  APF  iteratively  until 
the  total  loss  is  less  equal  to  Pg  within  tolerance.  This  is  the  optimal  solution. 

Step  -  8.  Determine  the  reduction  in  demand. 

LP  =  loss  of  step-1  -  loss  of  step-7/6. 

Step  -  9.. Determine  the  energy  loss  and  reduction  in  energy  loss. 

energy  loss  (E)  =  peak  power  loss  x  T  x  (Lf)^ 

LE  =  E  -  energy  loss  or  step-7/6. 

Step  -  10.  Determine  the  saving  due  to  loss  reduction  using  equation  7. 

b.)  Option-2:  Optimum  Location 
Step-1.  Same  as  Step-1  of  Option-1. 

Step-2.  Same  as  Step-2  of  Option- 1 . 

Step-3.  Determine  the  location  for  the  predetermined  capacitor  sizes.  This  is  done  as  follows; 
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3a.  Determine  the  normalized  reactive  current  fhi  for  the  ith  capacitor  using 
Ici=(  1 /a)  *  (Kvari/Kvart) 

where,  Kvari  is  the  ratings  of  capacitor  i  and  Kvart  is  the  base  kvar. 
i-1  0 

Icx=ZIcj,  ZIcj=  0 
j=l  j=l 

fhi=(Ici/2)+Icx 


(10) 


(11) 

(12) 


3b.  Compare  fhi  with  the  normalized  current(ficni)  of  Step-2.  The  location  for  the  ith  capacitor  is 
found  by  adding  all  the  section  lengths  upto  the  point  where  fhi=fxni. 


Step-4.  Call  Subprogram- APF  and  steps  8-10  to  determine  losses,  and  savings. 

One  of  the  requirement  of  this  algorithm  is  that  the  predetermined  capacitors  must  be  of 
different  sizes  successively  increasing  towards  the  source.  The  locations  determined  by  the 
algorithm  will  be  always  at  the  feed  points. 

c.)  Option-3:  Optimum  Sizes 

Step-1:  Same  as  Option-1. 

Step-2:  Same  as  Option- 1. 

Step-3:  Normalize  the  predetermined  locations  for  the  capacitors. 


Step-4:  Calculate  area  Ai,i+1  under  the  F(x)  curve  between  capacitor  locations  hi  and  hi-rl,  using 
hi 

Ai,i+1  =  J  F(x)  dx,  for  i=l,  2,...,n  (1^) 

hi+1 


Step-5:  Calculate  F(gi)  for  the  ith  capacitor  using 
F(gi)=(Ai,i+l)/(hi-hi+l) 

Step-6:  Calculate  the  size  of  the  ith  capacitor  using 

Qci=a*kvart*[F(gi)-F(gi-l)] 
and  F(go)=0 


(14) 


(15) 


Step-7:  Increment  i.  Return  to  Step-4  and  repeat  the  process  until  i=n,  where  n  is  the  number  of 
predetermined  locations. 


Step-8:  Call  Subprogram- APF  and  steps  8-10  to  determine  Losses  and  Savings. 
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4.  Test  Case 


Depicted  in  Figure  2  is  the  one  line  diagram  of  a  test  feeder  with  three  distinct  zones  of 
loading.  This  feeder  is  simulated  using  the  CapSize  computer  program,  with  a  substation  voltage 
of  12  load  factor  of  .95,  demand  charge  of  $160/lcw/yr,  energy  charge  of  15  mills/kwh  and 
capacitor  cost  of  $7/kvar.  The  feeder  data  file  is  in  Tables  1.  The  results  of  Capsize  simulation  are 
in  Tables  2  and  3  and  figure  3.  The  concentrated  load  results  refer  to  the  case  where  the  zone 
loads  are  concentrated  at  their  feed  points.  The  main  feeder  is  then  analyzed  as  one  radial  feeder 
with  no  zones.  This  case  gives  insight  to  into  the  performance  of  the  system  when  the  zone 
concept  is  dismissed. 


Table  1  -  Feeder  Data  File 


SN 

CONDUCTOR 

R 

X 

L 

p 

Q 

COMMENTS 

1 

477MCM/ACSR 

.216 

.4298 

.1477 

39 

29 

2 

477MCM/ACSR 

.216 

.4298 

.2424 

20 

15 

477MCM/ACSR 

.216 

.4298 

.1742 

13 

15 

4 

477MCM/ACSR 

,216 

.4298 

.0796 

13 

10 

5 

477MCM/ACSR 

.216 

.4298 

.1091 

82 

62 

6 

.216 

.1326 

53 

40 

7 

477MCM/ACSR 

.216 

.4298 

.1061 

62 

8 

477MCM/ACSR 

.216 

.4298 

.0379 

206 

9 

477MCM/ACSR 

.216 

'137 

103 

10 

477MCM/ACSR 

.216 

.4298 

.0758 

313 

11 

477MCM/ACSR 

.216 

.4298 

.0379 

39 

29 

12 

477MCM/ACSR 

.216 

.4298 

.0379 

11 

8 

0 

Z-2.  Lat-0 

2 

7 

0 

0 

2 

end  of  Zone-2 

1 

.216 

.4298 

.2803 

1276 

957 

2 

477MCM/ACSR 

.216 

.4298 

.0455 

49 

37 

3 

477MCM/ACSR 

.216 

.4298 

.1440 

33 

25 

4 

477MCM/ACSR 

.216 

.4298 

.0834 

116 

87 

5 

477MCM/ACSR 

WSMM 

.0834 

33 

25 

6 

477MCM/ACSR 

.2197 

67 

50 

l@l 

ISSWBWMIM 

0 

0 

3 

end  of  Zone-3 

477MCM/ACSR  .216 


477MCM/ACSR  .216 


477MCM/ACSR  .216 


477MCM/ACSR  .216 


477MCM/ACSR  .216 


477MCM/ACSR  .216 


477MCM/ACSR  .216 


477MCM/ACSR  .216 
477MCM/ACSR  .216 


Zone- 1  (Main  Feeder)  Q 


Legend  for  Table  1: 

SN  -  Section  Number 
R  -  Resistance  in  Ohms/Unit  Length 
X  -  Reactance  in  Ohms/Unit  Length 
L  -  Section  Length  in  Miles 
P  -  Real  Power  of  Load  in  KW 
Q  -  Reactive  Power  of  Load  in  KVAR 

Table  2  -  Capacitor  Size,  Location,  Loss  Reduction  and  Savings(Option  -  1) 


Zone  Number 


Concentiated  Load 


Zone  -  2 


Zone  -  3 


Zone 

Feeder) 


l(Main 


Distance  from 
Source  in  Miles 


2.3852 


3.3361 


2.66 

2.71 

3.02 


2.2716 


Size 

in  KVAR 


2861.7 


888.4 


77.5 

75.9 


Peak  Power 
Loss  in  KW 


83.4 


Peak  Power  Loss 
Reduction  in  KW 


Savings 
in  S/vT 


31.87 

8877.5 

1.23 

342.5 

0.65 

180.6 

32.21 

i 

8976.7 

Given  in  Table  2  is  the  theoretical  capacitor  sizes  and  locations  with  the  greatest  savings. 
Practically,  capacitor  banks  of  less  than  150  KVAR's  are  not  economically  feasible  to  install. 
Therefore,  some  of  the  smaller  banks  will  have  to  be  lumped.  Option  -  2.  is  used  for  this  purpose 
and  the  results  of  Option  -  2  CapSize  simulation  are  given  in  Table  3.  Figure  3  depicts  the 
improved  voltage  profile  due  to  installing  shunt  capacitors. 


4 

Table  3  -  Capacitor  Size,  Location,  Loss  Reduction  and  Savings(Option  -  2 


Zone  Number 

Distance  from 

Source  in  Miles 

User  Defined 
Size  in  KVAR 

Concentrated  Load 

2.3852 

2850 

Zone  -  2 

3.2982 

900 

Zone  -  3 

2.6655 

1050 

Zone  -  l(Main 

Feeder) 

2.2716 

1050 

Peak  Power 
Loss  in  KW 
83.4 


Peak  Power  Loss  Savings 
Reduction  in  KW  in  $/yr 

31.86 _ 8876.7 


1.22  341.8 


0.62 _ 174.3 


32.19  8970.5 
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Graph  Identification:  Starting  point 
Extreme  Left 

Section  Starts  at  Section  Ends  at 

1  9 

Comments 

Concentrated  Load 

After  9 

1 

12 

Zone-2 

After  12 

1 

6 

Zone  -  3 

After  6 

y 

1 

5.  Conclusion 

9 

Zone- 1  (Main  Feeder) 

.  f  M  ^  r _ 3 _ 

A  decorous  design  methodology  for  ohmic  loss  reduction  in  distnbution  feeders  is 
'presented.  Feeders  with  laterals  are  considered  through  the  zones  method.  A  computer  program 
called  Capsize  allows  three  design  options.  These  options  provide  a  powerful  tool  for  tailonng 
the  design  in  a  progressive  way  to  meet  practical  constraints. 
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